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Abstract-Support vector machine (SVM) is based on the special 
small samples theory with strong generalization ability, and is 
selected as an optimal theory for small samples classify problem. 
Electrical capacitance tomography (ECT) is a typical small 
samples and nonlinear mapping problem. In this paper the ECT 
image reconstruction algorithm based on SVM is proposed and a 
novel training method is proposed to improve the efficiency of 
SVM classifier by selecting active penalty parameters. 
Programming software based on Matlab6.5 and VC6.0, the 
simulation and experiment results indicates this algorithm has 
stronger space resolution and generalization ability, overall 
performance of the algorithm is better than some classic 
reconstruction algorithm such as LBP, etc. However, the 
algorithm requires the higher quality of training samples, 
otherwise less real-time, loss of application value. 
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Identification of Multiphase Flow widely in petroleum, 
chemical industry, river training and other areas of national 
product, as the complexity of multiphase fluid, random much, 
it is difficult using traditional detection techniques to 
accurately measure. Electrical Capacitance Tomography (ECT) 
technology is the use of multiphase medium tending to have a 
different dielectric constant, through the array of electrode 
capacitance changes, reflecting the multiphase medium 
pipeline distribution to extract the phase of each phase 
medium concentration, relative size, flow patterns and other 
characteristic parameters. The technology can be used to 
detect two-phase / multiphase flow processes, real-time 
display phase distribution images of multiphase fluid in the 
pipeline to achieve phase concentration, flow patterns and 
other characteristic parameters, is a measurement of high 
speed, non-invasive, low cost, applicable to a wide range of a 
new generation non-destructive functional imaging techniques [1]. 

Electrical capacitance tomography system includes 
capacitance sensor, data acquisition and signal processing unit, 
image reconstruction and analysis display three functional 
modules. As the "soft field" feature of system sensitive field, 



less independent capacitance measurements and other factors, 
makes the image reconstruction algorithm for the overall 
performance of the system plays a decisive role. As the 
current practical image reconstruction method are mostly 
qualitative algorithms (such as linear back projection 
algorithm), at the expense of accuracy to improve the 
real-time imaging, making the reliability of the extracted 
characteristic parameters less, and some high precision 
real-time imaging algorithm more poor, and no practical 
value. 

For less of independent capacitance value for the system, 
this paper design an image reconstruction algorithm for 
electrical capacitance tomography based on support vector 
machine (SVM) from statistical learning theory [2]. 

I. SVM Algorithm Design 

Support vector machine by Vapnik first proposed (Vapnik, 
1995), is the structural risk minimization method to achieve 
approximate. It integrates the largest interval of hyper-plane, 
Mercer core; convex quadratic planning, sparse solutions and 
slack variable technologies can be used for pattern recognition 
and nonlinear regression, which has become a standard tool 
for machine learning [3]. The 12-electrode capacitance sensor, 
for example, the system has 66 independent capacitance 
values, then an output unit for the support vector machines 
architecture shown in Figure 1, the number m of the input 
layer units are equal to 66, the number m, of hidden layer 
units are decided by support vector number in the system. The 
core of support vector machines makes the non-linear input 
space into linearly separable feature space through the inner 
core function K(x, x. ) [4] . 

Make {#>, (*)}"!, express a collection of non-linear 
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transformation from the input space to feature space, in which 
x is the input space vector. 



Bias b 




y„(d„) 



Output Layer 



Input Layer Hidden Layer 

Fig.l SVM architecture of output unit 

y k =d t =+1 Indicate existence of high dielectric 
constant material in the output unit. 

y t =d t =-1 Indicate existence of low dielectric 
constant material in the output unit. 

Where y k is output value of the output unit, d k is 
target value of the output unit. 

So the inner product kernel function as 

K(x, x, ) = (p(x) ■ p(x )) = JT Vj (x)<p l (x. ) ( 1 ) 

i=l,2,---,JV 
JV for training sample size , There is the number 
(x , d )f_ 1 of training samples. 

Selecting the radial basis function for the inner product 
kernel function : 

— — ii— — ii 2 

K(x,x.) = exp( — y x - x ) (2) 

Supposing the k output unit: 

A given training sample (x.,d.')", , Using convex 

quadratic planning algorithm for maximizing the objective 
function. 



MaxQ(a) = ^a i -i^^Va d^K^x.) (3) 

JV 

Constraints: ® / 'a,d, = 

® 0<«<C !=1,2,--,JV 
And to meet the constraints of the Lagrange 
multiplier cr , then its corresponding training sample {x } is 

a support vector. So that does not meet the constraints of the 
Lagrange multiplier set 0. C in the constraints © for the 
penalty parameter[5]. 



Corresponding to the optimal weight vector is 



W, = 



X a "' d ^^) 



(4) 



The output of the number k unit is as follows: 



y 



Za d K(x,x ) + b, 



(5) 



Where b t takes the b value arithmetic mean value of 
the support vector collection substitution equation (5). 

II. Image Reconstruction Simulation 

By the two dimensional finite element method analyzing 
system sensitive field, adopts triangle meshing and linear 
interpolation to seeking sensitive field electric potential 
distribution and examination pole plate capacitance [6]. If the 
sensitive area, the pipeline wall and the pipeline outer wall to 
the shielding layer region's meshing layer respectively is 4 
and 3 and 3. Programming software based on Matlab6.5 and 
VC6.0, meshing map shown in Figure 2. 
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Fig. 2 System meshing map and sample numbers 
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Simulation Capacitance value using equation (6) for 

normalized [7]. 

C" -C 

(6) 



C 



L +L 



-C , C 



L +L 



(7) 



" C'.-C 
If imaging medium is A,B and s A > s B ,among: C / , C'. is 

respectively the capacitance value when the tube filled with 



Where, L and L is the number of positive class and 
negative class point sets. C t and C are penalty parameters 
of positive class and negative class point sets. Then the 
equation (3) change as : 



No. 
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dynamic C 
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Fig.3 SVM image reconstruction comparison chart based on static and dynamic penalty parameter. 

A-phase and B-phase medium, C.™ is the capacitance value MaxQ(a) = V a. - \ V V a. a . d, d . K(x. , x . ) (8) 



when the tube filled with A, B for the two-phase medium. 

Based on 700 groups of training samples, select the 
equation (3) constraints © training network of the C = 8 to 
generate support vector set. Static C reconstructed image in 
figure 3 for the test samples (as distinct from the training 
sample) is the simulation images based on the support vector 
set through static threshold filtering. By the diagram shows, 
Space resolution, generalization ability of the algorithm is 
better than the linear back-projection, radial basis function 
neural network image reconstruction algorithm. But its space 
resolution error is still large, poor usability. 

As great differences of the numbers between positive 
class and negative class point sets. But the penalty parameter 
C in the equation (3) is static, this means that the number of 
certain types is more, then the penalty will be heavier, but we 
hope the same degree of punishment to positive class and 
negative class point sets, so structure the different penalty 
parameters : 



Constraints: ® Va,d. =0 



(2) 0<a, <C + , d, =1 

0<a<C_ , d. = -1 

i=l,2,---,JV 

Based on the same 700 groups training samples, select 
the equation (8) constraints (2) training network of the C = 8 
to generate support vector set based on the equation (7). 
Dynamic C reconstructed image in figure 3 for the test 
samples (as distinct from the training sample) is the 
comparison chart based on the support vector set through 
static threshold filtering. From the diagram shows, the space 
resolution error, generalization ability of the algorithm is 
better than image reconstruction algorithm based on static C , 
strong usability [8][9]„ 

III. Experiment and Analysis 
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Algorithm performance testing mainly study imaging 
experiments of the different numbers, different sizes and 
different locations of static objects, a single object, two 
objects, multiple objects, imaging experiments were selected 
to a diameter of 10mm and 20mm and 30mm of the Plexiglas 
rods. 



capacity, generalization ability, overall performance of the 
algorithm is better than LBP, etc. However, the algorithm 
requires the higher quality of training samples, otherwise less 
real-time, loss of application value. 



No. 



Type 



1 20mm 



2 20mm 



1 30mm 
1 20mm 
1 10mm 



Distribution map of the 
physical model 





Reconstruction images of 
the static C 






Reconstruction images of 
the dynamic C 






Fig.4 Imaging results test of plexiglass rods 



The following conclusions from the experimental results 
figure 4: 

CD The small diameter objects reconstructed images results 
rather vague, large diameter objects better reconstructed 
images; 

(D The closer to the center, the reconstruction image more 
blurred, the closer to the edge, the better the reconstructed 
image; 

(3) Images from SVM algorithm based on dynamic C 
than based on static C clearer, more precise. 

IV.CONCLUSIONS 

Image reconstruction problem of electrical capacitance 
tomography can be reduced to a class of linear inseparable 
problems, this paper based on dynamic penalty parameter, and 
convex quadratic planning in support vector machine, 
programming software based on Matlab6.5 and VC6.0, the 
image reconstruction of the simulation and experimental 
results show that the algorithm has stronger space resolution 
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